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ABSTRACT 

The experimental results of C0 2 adsorption and desorption in a packed column indicated that the concentration 
wave front at the center of the packed column differs from those which are close to the wall of column filled with 
adsorbent material even though the ratio of column diameter to the particle size is greater than 20, Mohamadinejad 
and Knox (2000). The comparison of the experimental results with one dimensional model of packed column shows 
that in order to simulate the average breakthrough in a packed column a two dimensional (radial and axial) model of 
packed column is needed, Mohamadinejad (1999). In this paper the mathematical model of a non-slip flow through 
a packed column with 2 inches in d.ameter and 18 inches in length filled with 5A zeolite pellets is presented. The 


comparison of experimental results of C0 2 adsorption and desorption for the mixed and central breakthrough of the 
packed column with numerical results is also presented 


INTRODUCTION 

The one-dimensional mathematical model approach to simulate a packed column indicates the lack of accuracy 
of concentration prediction away from the center of column, Mohamadinejad and Knox (2000). The experimental 
results on the concentration measurement show a substantial difference between the mixed and the central 
concentration in adsorption and desorption of C0 2 and H 2 0 on 5A Zeolite pellets. Numerous experimental data 
suggests that the channeling effect is significant even for a large ratio of column diameter to the particle size 
(d/d), Benenati and Brosilow (1962), Chu (1989), Vortmeyer and Michael (1985), Cohen (1981), and McGreavy 
(1986). Benenati and Brosilow (1962) indicated that channeling is very pronounce at the ratio of column 
diameter to the particle size of less than 20. Cohen (1981) showed that for the ratio of column diameter to the 
particle size of less than about 30 the a peak velocity of at about one particle diameter away from the wall; the 
velocity here ranged from 30 to 100% greater than the bulk velocity. The contribution to the radial temperature 
difference is mostly comes by different rate of adsorption along the cross section of the packed column which is 
based on the channeling effect in a nearly adiabatic packed column. In the adsorption of H 2 0 on 5A zeolite in a 
nearly adiabatic packed column the channeling effect is more pronounce for H 2 0 than CO, even after few hours 
when the cross-sectional temperature is constant, Mohamadinejad (1999). This observation is coherent with the 
fact that the heat transfer front is ahead of mass transfer front in a packed column. Therefore, a non-Darcian flow 
model (two-dimensional flow) has been developed to simulate the adsorption and desorption processes in a packed 
column. Finite-differencing numerical technique was used to solve the system of partial differential equations. For a 
packed bed the porosity varies with distance from the wall. Near the wall the porosity is higher than the bulk of the 
bed. This increases the permeability. A few particles away from the wall, the porosity equals the free stream value, 
Benenati (1962) and Roblee (1958). As a consequence of the porosity increase in the vicinity of the wall, the 
velocity of the flow parallel to the wall increases as the wall is approached and goes through a maximum before it 
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decreases to zero (to satisfy the no-slip condition). In general this leads to a net increase in flux, i.e., to the 
phenomenon called channeling, Nield and Bejan (1992). In this study, channeling effect on momentum, energy, and 
material balances was considered to be important enough so that two-dimensional adsorption in the packed bed 
must be modeled. 

Mathematical Model for Non-isothermal Multi-component Adsorption in a Packed Bed 

Momentum, heat, and mass balance equations can model the two-dimensional dynamic bed behavior. The 
mathematical model will be used to estimate the breakthrough curve for a certain constituent in the bulk gas. In 
return, this enables one to obtain the necessary parameters for predicting the transient behavior of the temperature 
profile and concentration of the gas for different initial parameters such as inlet concentration, temperature, and the 
fluid velocity. 

These equations were solved numerically by finite difference methods, namely the Newman methods (1968). A 
FORTRAN code was written to find the numerical solutions to the transient equations. 


Two Dimensional Adsorption Mathematical Model 

The complicated molecular diffusion of a component in a mixture is described by the Stefan -Maxwell 
equation. For the single component diffusion in a mixture, however, the diffusion coefficient D m j for the 
component is approximately related to the binary coefficients by the following relationship (Bird et al„ 1960) 


D m : = 


1 -y, 

X— 

^ n 


j*i 


( 1 ) 


For binary mixtures at low pressure, D[j can be estimated as it suggested by Slattery and Bird (1958). 
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Diffusion Model for Zeolite 


The rate of adsorption into the adsorbent pellets assumed to be approximated by the linear driving force 
approximation model, 


^r = k ef a s(<li -?i) 
ot 


( 2 ) 


Where k e f may be obtained by experimental procedure and a, is the interfacial surface area. The 
justification of assuming a linear diving force to model the adsorbed concentration in the solid phase has 
been well established by other researchers such as, Ruthven (1984), Do (1989), Grag and Ruthven (1972), 
and Sergent and Whitford (197 1 ) to name a few. 

Component concentration can be modeled as: 


dC, n d 2 C i d(uC t ) ' n d t rdC i ' (7-e;ag, 


dt ‘ ff “* dx 2 

Boundary /initial condition 


dx 


rdr dr 


dt 


fort<Oq=C i0 forO<x<L andO<r<R i 
fort>0 C] =q rXO at x=0 andO<r<R i 

fort>0 — =0 at x-L andO<r<R i 
dx 

()C 

fort =0atr=0andr = R 


The energy balance for the gas phase can be modeled as: 


(3) 


P*C 


g PS 


dT * h s a s (T g -T s ) (4) 


dt 


dx 2 


dx 


r dr dr 


Boundary /initial condition 
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at t<0, T g = T g0 for 0<x<L and 0<r<R i 
at t>0, T g = T 0x for x = 0 and 0<r<R t 
dT 

at t>0, — — = 0 for x - L and 0<r <R t 
dx 

dr 

at t>0, — — = 0 at r -0 for 0<x<L 
dr 


dT 

at t>0, k fr = h w (T w -T g ) at r=R forO<x<L 


The energy balance for the solid phase can be modeled as: 




(5) 


Boundary /initial condition 

at r< 0, T s =T s0 for 0<x<L and 0<r<R t 
at t>0, T s = T 0x for x = 0 and 0<r<R. 

at t > 0, — ^ = 0 for x-L and 0<r<R i 

dx 

dT 

at t >0, — - = 0 at r = 0 and r=R for 0 < x< L 

dr 

Rate of adsorption dqjdt can be substituted in the above equation. 

The energy balance for the wall can be written as: 
r)T 

PjC^^uRMT.-r.)- 2 *r,k<t, - T. ) 

Initial Condition 
at t<0, T w = T k/j 


(6) 


Bed Energy Equation Based On Effective Conductivity 

dT 

(£p g C pK +(l-£)p s C l)S )— = £k, 


d 2 T Kff.r a , rdT 


+ £ 




9/ " ’ dx’ r dr dr 




( 7 ) 


Boundary /initial condition 
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at t <0, T = T 0 for 0< x< L and 0<r<R i 
at t>0,T = T 0x for x = 0 and 0<r<R, 

at t >0, — = 0 for x = L and 0 <r < R, 
dx 

at t>0, — = 0 at r = 0 for 0< x<L 

dr 


at t>0, , k 


‘ITS 


ar 

dr 


= h.<T.-T\„ R ) M r = R, 


The governing momentum equation for cylindrical beds for fully developed flow is, from Vafai and Tien 


(1982). 


*.-p,Cu'-2-u+ 

dx W K r e dr dr 

Boundary Conditions 

K dP 


u = — 


/ u 3* 


at r-0 and u-0 at r= R i 


( 8 ) 


Where eis the porosity, and K and C are the permeability and inertial coefficient which are related on the 
porosity and the type of porous materials. In the above equation, the second term is the inertial effect which 
accounts for additional pressure drop resulting from interpore mixing found at higher Reynolds numbers, Valai 
(1982) and Ergun (1952). The third term is the Darcian force representing the pressure loss due to the presence of 
solid particles. The last term is the viscous shear force representing the resistance to the flow caused by sheer stress 
along the solid boundary. This term accounts for the no-slip boundary condition at the solid boundary. In this study 
the entrance effect is not considered since the flow is fully developed after one to two-particle distance from the 
entrance, see Vafai and Tien (1982). 


Porosity Variation 

£ = e„\l + aexp(-by/d p \ (9) 

Where £00 is the free-stream porosity, y is the distance from the wall, d is the particle diameter, and a is taken to 
be 1.4, Nield and Bejan (1992). b is experimental parameters that depend on packing and particle size, it varies from 
2 to 8. 
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The empirical coefficients K and C, which are given by the relations developed by Ergun (1952) for flow in a 
packed bed: 


a: = 


d 2 e 3 


150(l-e) 2 


C = 


1.75(1 -e) 


( 10 ) 


(ID 


de- 


The variable C and K, both are function of the bed porosity and particle diameter, d. The porosity in a packed 
bed increases from the center of bed, free-stream porosity, to a maximum of one at the bed-wall boundary. This 
increase is confined within few particle diameters from the wall, Benenati (1962), Roblee (1958). 

In the above two dimensional equations, the term which represents the radial diffusion, is 


1 3 frdC' ■ 
r dr ( dr . 

Where C is a variable. By carrying out the derivative, it can be recast into 

\_dc_ tfc_ 

■r dr dr 2 

At the center where r = 0, the first term is not finite. But 
/ 18C > \ d 2 C 


( 12 ) 


(13) 


Limit 


r- ► 0 


r dr 


dr 2 


(14) 


By L'Hospital’s rule. Therefore the term, eq 12, in 2-dimensional form for the center point is replaced by 


IJL 

r dr 


< 'rdC ^ 

dr 


= 2 


d 2 C 

dr 2 


(15) 


Therefore, the diffusional term in the discretized forms of two-dimensional PDFs at center grid is replaced by 
eq (15). 


Calculation of Thermal conductivity for 2.D Flow 


In this study two different equations were used to calculate the effective conductivity in the packed bed. One is 


7 



based on the works of Kunii and Smith (1960), and the other one is based on the experimental work of Fahien 
(1954). 

The effective thermal conductivity in the axial and radial direction, k e ff iX an d k ejf„r are related as, 

k 'ff,x =k ° +k f.x ( 16 ) 

k eff.r = k ° + k f,r 

In the above equations, the radial and axial conductivity is the combination of two terms. The First term is the 
stagnation conductivity, while varies from a bulk conductivity to fluid conductivity with distance from the center to 
the column wall. Therefore, it depends on the porosity variation, which also is a function of bed parameters. The 
second term is due to the dynamic or dispersion conductivity, which incorporates the mixing, causes by flow 
through the particles. This conductivity can be calculated by theoretical equations. Based on the work of Kunii and 
Smith the following equation result. 


Stagnant Conductivity k„ 

Kunii and Smith (1960) presented theoretical equations for estimating the stagnant conductivity, k 0 . The 
stagnant conductivity can be found If ( 5 , solid thermal conductivity, is given. 

— = e + (\-£)/^> + (2/3)(k f /k s )] < 17 > 

k f 

Where kf is the thermal conductivity of fluid and <f> is the contribution of solid to solid heat transfer through 

fluid film around a contacting point of neighboring particles. 0 Is given by 

<p=<p 2 +(<l> l -0 2 )[(e ~.26)/.216] 

for .476 >£> .26 

0=01 (18) 

for e > .476 

0=0 2 
for £ < .26 

Where 0, and <p 2 are given in a schematic form, and are being interpolated linearly in tabular form in the 
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computer program. 


Effective Radial Conductivity kf >r 


The thermal conductivity in radial direction for packed bed is given by Baron (1952) as 

C„P,^k,^N r , t =8- , 0 (19) 

Where N FtH is Peclet number, therefore the effective thermal conductivity in radial direction would be Yagi and 
Kunii (1957). 


= •— + (a/3 )N Rep /V Pr 

k f k f 

where ( ap ) = 1/N Pe = . 1 to . 125 


(20) 


Effective Axial Conductivity keffa 


A similar equation can be derived for effective thermal conductivity in the axial direction, Yagi, Kunii, and 
Wakao, (1960). 

w.* 


k,ff r ko 

- — H A/V n . N p. 

I L Ke f n 

Kf Kf 


( 21 ) 


where A = .5 to 1 .0 

Incorporating the effects of porosity variation into the effective conductivity, the effective conductivity reduces 
to. Hunt (1987) 


= (l + a~ exp(-bR))^— + (aP)^^-N Re Np t 

k f k f d 

where ( aP )= 1/N PeH = .1 to .15 


( 22 ) 


And in the axial direction 

k 


7 ^-=(l + a- cxp(-bR))^~- + X^^-N Rep Np T 


k f d 


(23) 


where A = .5 to 1.0 

Where a is chosen such that kf/kf equals one at the wall as velocity becomes zero. The variation of dispersion, 
mixing length, is 
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^ - R t - r for R t - r<] 
d 


( 24 ) 


Where Rt is tube diameter. An expression similar to porosity variation was used to predict the mixing length 
variation. Hunt (1987). 


Calculation of Mass diffusivity For 2.D Flow 


A similar theoretical approach can be taken for the calculation of diffusivity in the radial and axial direction. 


Effective Diffusivity 

Effective diffusivity follows the same expression as in thermal conductivity 

D eff , x = D. + D fx 

D eff,r = D ° +D f,r 


(25) 


Effective Radial Diffusivity Df >r 

The effective diffusivity in the radial direction by analogy to heat transfer is 

£g± = J£- + ( a p)N K ' F N sc , 

Df U f 

where ( aj3 ) = 1/N PeH — .1 to . 125. 

Effective Axial diffusivity D e ff yX 

A similar equation can be derived for effective diffusivity in the axial direction, 


( 26 ) 


D eff.x D 


D f D f 


= ~z t" N SC i 


(27) 


where X = .5 to 1 .0. 

Incorporating the effects of porosity variation into the effective diffusivity, the effective diffusivity reduces to, 
Hunt (1987) 
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( 28 ) 


1 = (1 + a' exp (~bR))~+ (af5)^-N R N^, 
D f D f a 

where ( a/3 ) = 1/N Pefi - A to . 15. 

And in the axial direction, 

= 0 + a ' exp (-bR)) ~ + A N Re/> N Sc , 

L) j U y d 

where A = .5 to 1 .0. 


Effective Radial Thermal Conductivity Based on Fahin Equations 


In contrast to weak effect of mass diffusion on radial mass fraction, the thermal conductivity profile has a 
strong effect on both temperature and mass adsorption. A thermal conductivity profile for different ratio of dp/d t 
was approximated by Fahien (1954), 


k eff.r 


k 


M ~ 


2 3 

=< k > (k* + 3^ -KY-J+ 2 (k m M - k’ 0 )-Zj) 0 < r <r M 


M 


' M 


-< k > {k M - ( k M - k w )- — ) r > r } 

1 — 


M 


3-0 -.9 k 0 r M -k w 


2 r M 3 - 5r M 2 + 2 


1 -n 


M 


1 + ^ +. 1 ^' 


(30) 



a 

a = d t I d p 

Where K / is the effective radial conductivity at the center of the column wall, K M * is the maximum effective 
thermal conductivity, K w * the effective thermal conductivity near the wall, <tf>is the average effective thermal 
conductivity, and r is the location of maximum in conductivity profile. These conductivities are obtained the Argo 
and Smith (1953) equation using the void fraction values as a function of radial position. According to Argo and 
Smith 


k €ff -£\ 


, d p C pg G ( 
k g + — + 4 


N p e E 


— V(0.173)(r a 3 /100 4 ) 


2-a 


(l-e) 


hk s d p 


2 k s +hd 


(31) 


Where O is the emissivity of solid particle, T a j s the average temperature. In the above equation 
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h — h c + h r 4- hp 

K =\.95C pg GN Pr - 2n N Re ~ 5 ' N Re < 350 
K = 1.06C w GyV Pr - 2/3 /V Re - 41 N Re > 350 

k r (2k s +hd ) 

«r = 
k r = 


log, 0 k p = 1.76 + .0129^ s /a 

Effective Radial diffusivity Based on Fahin Equations 


d p k s 




kp{ 2*, +wj 


2 3 

Deff.r =<D> (D* + 3(d; -D*)-^ t + 2(d^ - D 0 * J-^) 0 < r < r M 


* M 


r M 


Dfff,r —<D> (D 


M 


'M 


1 -rj 


) r > 


M 


(32) 


Af 


Oa, = 


3.o -.9 p 0 y 

l+'W +-1'M 2 


D w * is obtained from work of Fahin and Smith (1955) to be 
D’ a =9/8V 0 (l + 4.85cT 2 ) 


Where V is velocity at the center of packed bed. 

The effective thermal conductivity in the wall layer of thickness Rp-dp/2 , k ew , is defined and h w is considered 
as a correction factor based on the difference k e >ff\ r and k €W (Kunii and Suzuki, 1966). 

1 


b k 

K ew 


V 


1 


1 


(33) 


a H< Npr/V R e h^dplkj 

Where Of w denotes the contribution of fluid mixing in the wall layer and is taken as 0.2. h w represents the heat 

transfer coefficient of the thermal boundary layer, which develops on the wall surface. This becomes dominant at 
high Np e p and is given by Blasius type equation as (Suszuki, 1990), 
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h' w d p !k f =c(N Pl U3 N Rep V4 ) 


(34) 


Where C is an experimental coefficient with value of 0.1 to 0.2 (Kunii, Suzuki and Ono, 1968). k C wo' s obtained 


by the following equation similarly to eq 17. 



£„+(! — ej/k +(2/3X*/ /*,)] 


(35) 


Where £ w denotes void fraction in the wall layer of about 0.7. 


NUMERICAL SOLUTION 

The solution to the non-equilibrium , non-isothermal adsorption/desorption problem must be done 
numerically. In this study , for a two component mixture , the numerical model would involve the solution of 
several coupled differential equations: one mass balance equation , two mass balance of rate equations (solid 
phase), one total mass balance equation , one momentum equation, one heat balance of fluid flow equation, one 
heat balance of solid equation , and one equation for heat balance for a wall 

The finite difference technique , which is mostly used for processes with varying boundary and initial 
conditions, is a more convenient method to use for solving the coupled, partial differential equations (PDE). In 
this work , the PDE's were discretized by first or second order differences in time and spatial dimensions. The set 
of discretized finite difference equations was solved simultaneously by the implicit method. Based on the stiffness 
and the sharpness of momentum equations and the mass and heat transfer fronts, the implicit method of 
Newman (1968,1967) was used to simulate the adsorption process. 

The examination of convergence or the rate of convergence of the numerical methods used in this study 
relied on actual testing of the iterative methods. Since the set of equations is both coupled and nonlinear, an 
analytical procedure for determination of convergence can not be used . 

In Newman *s method, the number of grids could be reduced to 50 grids and the time step ranges adjusted 
from 12 to 60 seconds. The iteration was stopped if ( C m+t -CJIC m +, was less than 1.0E-4 for each grid point. 

The stability and the rate of convergence in these PDE's are mostly affected by rate of adsorption, which are 
related to isotherms. At low concentrations , the adsorption affinity of N 2 on 5A Zeolite is less than C0 2 by about 
an order of magnitude. The adsorption affinity of CQ 2 is also less than H 2 0 by about an order of magnitude , 
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Mohamadinejad (1999). The steep concentration gradient of HjO relative to C0 2 and C0 2 relative to N 2 causes 
the rate of convergence to be smaller for H 2 0 in comparison to C0 2 and smaller for C0 2 in comparison to N r For 
single component adsorption of H 2 0 and C0 2 on 5 A Zeolite , the time step for an HjD run must be much less than 
the time step used for a C0 2 run , otherwise the numerical model would not converge , i.e.. Also, the number of 
iterations for convergence with a given time step is much larger for Hft than CO p Mohamadinejad (1999). The 
second most important criteria is the inclusion of a diffusion rate into the mass and energy balance equations. 
The implementation of diffusion coefficients makes the PDE's parabolic functions and in turn causes the 
equations to be more stable, this is because diffusion exercises a smoothing effect on the PDE’s, (Sewell, 1988). 


Test Apparatus 

Verification of the numerical models is accomplished with test data from the Molecular Sieve Material Bench 
Test (MSMBT), a Marshall Space Flight Center apparatus. The test apparatus consists of a small packed column 2 
inches in diameter and 10 or 20 inches long, depending on the configuration. Instrumentation of the column 
includes temperature probes and sampling tubes for measurement at sorbent material endpoints and intermediate 
points. Continuous measurements of the exit CO 2 partial pressures, and all temperatures are recorded. A gas 
chromatograph is used to capture the breakthrough at interior bed locations. The MSMBT is fully described 
elsewhere ( Mohamadinejad, 1999). Any significant changes made in the test configuration will be discussed in the 
following sections as appropriate. 

The MSMBT is used to obtain data empirically that is not available otherwise or not reliable in the open 
literature, such as heat of adsorption and lumped mass transfer coefficients. Testing to obtain empirical data is 
performed in such a way to isolate the phenomenon of interest as much as possible. Heat transfer coefficients, (not 
available otherwise due to the use of unique NASA foam insulation for superior adiabatic conditions) for example, 
were obtained by flowing dry, heated nitrogen through a sorbent bed desorbed in the presence of nitrogen. Nitrogen 
gas was used since it is the main carrier gas. The adsorption affinity of nitrogen on 5A zeolite is negligible at 
high temperature, therefore the heat generated during the dry-heated nitrogen through the column is minimal. 
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Also, the curve fitting of experimental and model results was done at the later portion of curve where the 
temperature is high . However, the obtained heat transfer coefficients fits the entire curve fairly accurate , 
Mohamadinejad (1999). 

The second use of the MSMBT is to verify, following ascertain of empirical values, the accuracy of the model 
under conditions similar to those expected in the International Space Station Alpha (ISSA) . This consists of 
comparing a series of MSMBT runs at varying conditions to the results of the computer model ran at the same 
conditions. Of special interest is the capability of the model to predict transient bed temperature and partial pressure 
for conditions not used to obtain the empirical values. 

COMPUTER MODEL VERIFICATION 
Carbon Dioxide Single Component Adsorption 

The test results of will be used to verify the two dimensional bed models for carbon dioxide adsorption. Test 
specifications are shown in table 1 . 

Two phenomenon were observed to occur during this adsorption test, which are not generally accounted for in 
mathematical models of this type. 

The first phenomenon is the significance of bed channeling in the 1.875 inch ID cylindrical column with 
sorbent pellets of 8-12 mesh (app. 2 mm in diameter). The existence of significant bed channeling is evident in 
Figure 1, which shows the breakthrough of CCb in the column. Note the discrepancy in the breakthrough at the 
material exit centerline (labeled "Outlet") from a point 5 inches downstream of the material exit ("Outlet Beads"). 
The downstream point measures all gas exiting the column mixed via turbulence through 5 inches of 3-mm glass 
beads. Consistent with the assumption of flow channeling at the column wall, the breakthrough is earlier for the 
mixed gas than that at the centerline of the column. The two-dimensional model was developed to account for the 
channeling, and to aid in derivation of a technique to intelligently apply a channeling factor to the more CPU 
efficient single dimension models. 

The second phenomenon observed is the importance of including nitrogen co-adsorption for accurate modeling 
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of the carbon dioxide breakthrough with nitrogen as the carrier gas. As it was discussed in previous paper, 
Mohamadinejad and Knox (2000). 

Model inputs are shown in table 1. Pellet Specific Heat, is based on Davison Chemical Co. supply. 


Table 1 . Model Inputs for Two Dimensional Flow 


Flow rate, 

28.04 standard, liter/min 

Inlet Pressure, 

15.587 Psia 

C0 2 Partial Pressure, 

6.14 mm Hg 

H 2 0 Partial Pressure, 

6.0352 mm Hg 

Initial Temperature, 

77 T 

Bed Void Fraction, 

.35 

Interfacial surface 

area, 635 

Pellet Density, 

75. Lb/ft 1 

C0 2 Mass Transfer Coef, 

.02 ft/hr 

C0 2 Heat of Adsorption, 

18000.0 Btu/lb-mole 

N 2 Heat of Adsorption, 

8988 Btu/lb-mole 

Specific Heat (Pellet), 

.2 Btu/lb-°F 

Thermal conductivity (pellet), 

. 1 Btu/ft 2 -hr-°F 


Two Dimensional Model Verification 


The determination of input values, mass transfer coefficients and verification of the two-dimensional single and 
multi-components material flow-through adsorption and desorption model, is described in this section. 

Carbon Dioxide Adsorption 


Results of the model comparison of C0 2 single component are shown in Figure 2 and 3. As is evident, the 
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comparison is very favorable. Note that two sets of model data are presented; averaged data and centerline data. 
The centerline data is the central node radially and the last bed material node axially. The simulated two- 
dimensional centerline data matches the one-dimensional simulated results , Mohamadinejad and Knox (2000). 
Averaged data is also at the last bed node axially, but is an average of all the radial nodes. The averaged data is thus 
representative of gas after mixing in the glass beads, or the test data labeled "Outlet Beads". The centerline data is 
representative of test data taken at the center of the bed radially and at the end of the bed material axially. 

Based on the above discussion, it is expected that centerline data - the line with filled markers in Figure 1 will 
compare with the small circular markers. As seen from the figure, this comparison is indeed favorable. Average 
data should be compared with the small triangular markers, and once again, this is a favorable comparison. As a 
result, it can be inferred that the two-dimensional model does correctly model the channeling observed in this test. 

As seen from Figure 3, the temperature comparison is also favorable. The increase in fidelity of the two- 
dimensional model is evident by comparison of Figure 4 with the results from the one-dimensional model, 
Mohamadinejad, Knox (2000). The two-dimensional simulation more closely follows the actual temperature peak, 
both in time and in magnitude. 

The short time of breakthrough allows it to use Ideal Adsorption Solution Theory (IAST), even though it takes 
more CPU time than using the Langmuir isotherm. A mass transfer coefficient of .017 ft/hr was used in predicting 
the CO 2 breakthrough curve. 


Carbon Dioxide Desorption 

Results of the model comparison are shown by the solid lines in Figures 5 and 6. The model prediction of 
centerline and average breakthrough is definitely matches the obtained experimental data. The temperature profile 
result of the two-dimensional model also estimates the experimental data fairly well. The few degree discrepancy 
between the model and experimental data is largely due to predicting the heat transfer coefficients between the 
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packed bed and the wall and between the wall and the surrounding. The two parameters have strong effect on radial 
temperature profile. Trial and error must be performed in order to obtained the corrected values. The IAST was used 
to predict the CO 2 /N 2 mixture isotherm. A mass transfer coefficient of .017 ft/hr was used for the prediction of 

breakthrough curve. 


COjTemperature and Concentration of Model Results in Radial Direction 

The influence of porosity variation and the no-slip flow on temperature and concentration fronts of C0 2 
adsorption along the radial direction is significant as it shown in Figure 6 and 7. The early breakthrough of C0 2 
increases from the central line up to the wall of the packed bed column. It is evident that the dispersion in axial 
direction is lesser toward the wall than the center of the column, specially at beginning of the breakthrough. Figure 
8 and 9 show the model results of temperature variation at outlet of adsorbents in the radial direction. Even though, 
the effort made to make the column adiabatic the temperature variation in the radial direction is significant. 
Although, the larger porosity near the wall causes an early breakthrough of the C0 2 , but the lower temperature at the 
wall makes the C0 2 to be adsorbed more because of higher affinity of adsorption at lower temperature. It should be 
noted that the porosity variation is a decaying exponential phenomenon and therefore more pronounce within 10 
percent of the wall radius. However it occupies where the surface fraction is the largest. 

Figure 10 and 11 presents the concentration breakthrough for C0 2 desorption for different points along the 
radial direction. The effect of porosity variation and the no-slip flow on stripping the CO : is that it takes longer to 
strip the center of the column from C0 2 than close to the wall. Figure 12-13 shows the model results of 
concentration and temperature variations at outlet of the sorbent with respect to radial position. It is also evident that 
because the centerline is more condense the temperature is lower because of the effect of heat of desorption. 

CONCLUSIONS 

Based on analytical and experimental investigation of two-dimensional connective flows in porous media the 
following conclusions are drawn: 

• The experimental results from the laboratory scale-fixed bed adsorber are quantitatively consistent with the 
one-dimensional model at the column center. The average concentrations of cross sectional bed obtained by test 
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result deviates from the column center concentration in one-dimensional model appreciably and are fairly good with 
two-dimensional model. This indicates the strong effects of porosity variation along the radial direction of column 
bed on the temperature, concentration, and velocity field can be modeled by none-slip flow model. 

• A linear driving force mass transfer model for the solid phase provides a reasonable fit to experimental 
adsorption data. 

• The average and centerline C0 2 concentration breakthrough and temperature profiles of the test results matches 
the two-dimensional model reasonably well. 

• It is evident that in modeling the packed bed column with one dimensional flow, an average velocity and 
porosity should be considered so that at least some of the effects of radial variation of these two parameters are 
introduced into the one dimensional model. 


NOMENCLATURE 

2 3 

a Surface area of pellets per unit volume of pellet ft /ft 

th 3 

C Gas phase concentration of i component in the bulk, lb mole/ft 

C Constant in Darcy equation 

, 111 

C jo Gas phase concentration of i component at boundary or initial condition, lb mole/ft 3 
C s Heat capacity of gas phase, Btu/lbm-R 
C ps Heat capacity of solid particle, Btu/lbm-R 

C pn Heat capacity of column wall , Btu/lbm-R 

d, d' Column diameter, ft 

d p Particle diameter, ft 

D Diffusivity, ft /hr 

D o Stagnent diffusivity, ft /hr 

£> Axial diffusion, ft /hr 

D Binaiy Molecular diffusion , ft /hr 
D Molecular diffusion in mixture, ft /hr 

mi 

2 

D f Fluid flow diffusivity, ft /sec 
D Effective axial diffusivity , ft /sec 

D Effective radial diffusivity, ft /sec 

D, Axial diffusivity of fluid flow, ft /sec 

D fr Radial diffusivity of fluid flow, ft /sec 


19 



d Bed diameter, ft 

2 

G Superficial mass velocity, lb m /ft hr 

2 

h Effective heat transfer coefficient for column insulation, Btu/ft -hr 
° 2 
h Heat transfer coefficient between the gas and the column wall, Btu/ft -hr 
■* 2 
h Heat transfer coefficient between the gas stream and the sorbent , Btu/ft -hr 

K Constant in Darcy equation 

k f Fluid flow axial conductivity, dispersion, Btu/ft-hr-R 

k Stagnant thermal conductivity, Btu/ft-hr-R 

0 

k Effective axial conductivity, Btu/ft-hr-R 

efix _ 2 

k a Effective overall mass transfer coefficient, monodisperse particle, ft /hr 
k g Effective radial conductivity, Btu/ft-hr-R 

k f Axial conductivity of fluid flow, Btu/ft-hr-R 
k f Radial conductivity of fluid flow, Btu/ft-hr-R 
k Solid thermal conductivity in axial direction, Btu/ft-hr-R 

s.x 

k Solid thermal conductivity in radial direction, Btu/ft-hr-R 

s,x 

Jt_ Fluid film mass transfer coefficient of i™ component, ft/hr 

k Effective mass transfer coefficient of i™ component, ft/hr 
L Column length, ft 

M Molecular weight of adsorbate i, lb/lb mole 
n Number of component 

N p Peclet number, cpgPgu/k 

N Reynolds number, dpPgu/p 

N Schmidt number, p/Dpg 

2 

P Total pressure, mmHg or lbf/ft 

2 

P Partial pressure of component i, mmHg or lbf/ft 

3 

q Amount adsorbed in the solid, lb moles/ft of solid 

* lh 3 

q, Solid phase concentration of i component in equilibrium with gas phase, lbs moles/ft of 
solid 

3 

q i Volume average solid phase concentration of component i, lb moles/ft of solid 

r Radial position, ft 

3 

R Ideal gas constant 555 mm Hg ft /lb mol R 
R Inside wall diameter of column, ft 

1 

R Outside wall diameter of column, ft 

o 

t Time , hr 

T Temperature R 

T Ambient temperature, R 

o 

T gas temperature, R 

T Wall temperature, R 
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T Solid temperature, R 

u Interstitial velocity, ft/hr 

x Axial position 

th 

x Mole fraction of i component in the solid phase 

' th 

y Mole fraction of i component in the gas phase 

■f 

GREEK 

a/J Constant in effective conductivity equation 

e External bed void volume 

X Constant in effective conductivity equation 

<j) Constant in stagnant conductivity equation 

3 

p Density of gas phase, lb mole/ft 

PS 3 

p Density of solid phase, lbs/ft 

p Density of column wall, lb mole/ft 

H Fluid viscosity 

a Emissivity 

AH Heat of adsorption, BTU/lb of solid 
SUBSCRIPTS 


e Effective 

eff Effective 

/ In the fluid phase 

lh 

i i component 

o Outside, initial 

pg Gas phase 

s Surface 

t Total 

w Wall 

SUPERSCRIPTS 


Average value 
Equilibrium value 
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Figure 1 . C0 2 Breakthrough for Various Radial Positions 
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Figure 2. CO, Breakthrough Comparison Test with Two-Dimensional Model Results 
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Figure 4. CO2 desorption Comparisons Test with two-dimensional results 
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Figure 5. C0 2 Gas Temperature Desorption Comparison Test with Two-Dimensional Results 
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Figure 6. Calculated CO2 Breakthrough at the Outlet of the Bed for Different Radial points 
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Figure 7. Calculated CO 2 Adsorption Concentration at the Outlet of the Bed for Different Time Along the Radial 
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Figure 8. Calculated Adsorption Temperature at the Outlet of the Bed at Different Radial Location 
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Figure 9. Calculated Assorption Temperature at the Outlet of the Bed for Different Time Along the Radial 

Direction 
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Figure 10. Calculated CO2 desorption Concentration at the Outlet of the Bed for Different Radial points 
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Figure 11. Calculated CO 2 Desorption concentration at the Outlet of the Bed for Different Time Along the Radial 

Direction 
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Figure 12 . Calculated Desorption Temperature at the Outlet of the Bed at Different Radial Location 
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